Astronomy & Astrophysics manuscript no. RRLyr'par' Kolenberg astroph 


© ESO 2010 


April 30, 2010 





An In-Depth Spectroscopic Analysis of the Blazhko Star RR Lyr * 
I. Characterisation of the star: abundance analysis and fundamental parameters 

K. Kolenberg 1 , L. Fossati 2 , D. Shulyak 3 , H. Pikall 1 , T.G. Barnes 4 , O. Kochukhov 5 , and V. Tsymbal 6 

1 Institut fur Astronomie, Universitat Wien, Tiirkenschanzstrasse 17, 1 180 Wien, Austria, 
e-mail: katrien.kolenberg@univie . ac . at , holger .pikall@univie . ac . at 

2 Department of Physics and Astronomy, Open University, Walton Hall, Milton Keynes, MK7 6AA, UK. 
e-mail: 1 . f ossati@open .ac.uk 

3 Institute of Astrophysics, Georg- August-University, Friedrich-Hund-Platz 1, D-37077, Gottingen, Germany, 
e-mail: denis . shulyak@gmail . com 

4 The University of Texas at Austin, McDonald Observatory, 1 University Station, C1402, Austin, Texas, 78712-0259, USA. 
e-mail: tgb@astro . as . utexas . edu 

5 Department of Astronomy and Space Physics, Uppsala University, 751 20, Uppsala, Sweden, 
e-mail: oleg@astro .uu.se 

6 Tavrian National University, Vernadskiy's Avenue 4, Simferopol, Crimea, 95007, Ukraine, 
e-mail: vadim. tsymbal@gmail . com 



ABSTRACT 

Aims. The knowledge of accurate stellar parameters is a keystone in several fields of stellar astrophysics, such as asteroseismology 
and stellar evolution. Although the fundamental parameters can be derived both from spectroscopy and multicolour photometry, the 
results obtained are sometimes affected by systematic uncertainties. In this paper, we present a self-consistent spectral analysis of the 
pulsating star RR Lyr, which is the primary target for our study of the Blazhko effect. 

Methods. We used high-resolution and high signal-to-noise ratio spectra to carry out a consistent parameter determination and abun- 
dance analysis for RR Lyr. The LLmodels code was employed for model atmosphere calculations, while SYNTH3 and WIDTH9 
codes were used for line profile calculation and LTE abundance analysis. We provide a detailed description of the methodology 
adopted to derive the fundamental parameters and the abundances. Stellar pulsation reaches high amplitudes in RR Lyrae stars, and as 
a consequence the stellar parameters vary significantly over the pulsation cycle. The abundances of the star, however, are not expected 
to change. From a set of available high-resolution spectra of RR Lyr we selected the phase of maximum radius, at which the spectra 
are least disturbed by the pulsation. Using the abundances determined at this phase as a starting point, we expect to obtain a higher 
accuracy in the fundamental parameters determined at other phases. 

Results. The set of fundamental parameters obtained in this work fits the observed spectrum accurately. Through the abundance anal- 
ysis, we find clear indications for a depth-dependent microturbulent velocity, that we quantified. 

Conclusions. We confirm the importance of a consistent analysis of relevant spectroscopic features, application of advanced model 
atmospheres, and the use of up-to-date atomic line data for the determination of stellar parameters. These results are crucial for further 
studies, e.g., detailed theoretical modelling of the observed pulsations. 
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1. Introduction 

The modelling of pulsational signals requires the knowledge of 
stellar parameters and primarily accurate values of the effec- 
tive temperature (T e g) and metallicity (Z). The determination 
of fundamental parameters can be performed by different meth- 
ods (some examples for stars from B- to G- type are Fuhrmann 
et al. 1997, Przybilla et al. 2006, and Fossati et al. 2009) that 
do not always lead to consistent results. Thus, it is important to 
choose a methodology that allows us to constrain the parameters 
of the star from the available observables (usually photometry 
and spectroscopy) in the most robust and reliable way. 

RR Lyr is the prototype and eponym of its class of pulsat- 
ing stars. RR Lyrae stars play a crucial role as distance indi- 
cators. Their evolutionary stage (He burning in core, H burn- 
ing in shell) makes them useful tracers of galactic history. These 

* Data obtained with the 2.7-m telescope at McDonald Observatory, 
Texas, US 



classical pulsators display radial oscillations (the simplest type 
of pulsation) with large amplitudes which makes them useful 
touchstones for theoretical modelling. RR Lyr is one of the 
best studied stars of its class. Almost a century ago Shapley 
(1916) discovered that it shows a strong Blazhko effect, i.e. 
a (quasi-)periodic modulation of its light curve shape in am- 
plitude and phase. The Blazhko effect in RR Lyr has been 
closely followed over the past century, and changes have been 
reported both in the strength and the duration of its Blazhko 
cycle (Szeidl 1988, Kolenberg et al. 2006). Some well-studied 
stars even show multiple (variable) modulation periods (see, e.g., 
LaCluyze et al. 2004). Despite numerous attempts to model the 
phenomenon, the Blazhko effect has eluded a satisfactory ex- 
planation so far. Recently obtained high-precision photometry 
from ground-based or space-borne precise instruments also in- 
dicate that Blazhko modulation may be a much more common 
phenomenon than initially thought: as many as half of galactic 
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RRab stars may be modulated (Jurcsik et al. 2009; Szabo et al. 
2009; Kolenberg et al. 2010). 

In order to constrain the viable models for the Blazhko effect, 
it is vital to obtain accurate values for the fundamental parame- 
ters (and their variations) for modulated and non-modulated RR 
Lyrae stars. This has been the main motivation for the study pre- 
sented in this article. 

RR Lyr is the only star of its class to have a directly de- 
termined parallax, recently measured with the HST/FGS, by 
Benedict et al. (2002), to be Tr(FGS) = 3.82 + 0.2 mas (d = 262 
+ 13 pc). When a small ISM correction of A v = 0.07 is applied, 
this new distance results in an Mf R +0.61/2 'Jj mag which 
corresponds to ^ 49 ± 5L Q . 

Fundamental parameters of RR Lyr have been obtained by 
several authors with a variety of methods (e.g., Lambert et 
al. 1996; Manduca et al. 1981; Siegel 1982; for a summary 
see Kolenberg 2002). The published fundamental parameters 
of RR Lyr display a considerable range both in and log g 
due to the large pulsation amplitudes. According to these anal- 
yses, RR Lyr's T e g varies over its 13h36min pulsation cycle be- 
tween approximately 6250 and 8000 K and its logg between 
2.5 and 3.8 (extreme values). Superposed on the large variation, 
the Blazhko cycle leads to an additional variation of the fun- 
damental parameters. Jurcsik et al. (2008) recently showed that 
also the mean properties of modulated RR Lyrae stars change 
over the Blazhko cycle. Element abundances of RR Lyr were 
obtained previously by, e.g., Clementini et al. (1996), Lambert 
et al. (1996), and Takeda et al. (2006). 

The main goal of the present work is to perform a self- 
consistent atmospheric and abundance analysis of RR Lyr 
that reproduces all of its photometric and spectroscopic data. 
Furthermore, we want to investigate the degree to which the de- 
rived fundamental parameters depend on the applied methods. 
Considering the structure of the available models, especially the 
position of the convective zones and the zones of nuclear fusion, 
the measured abundances of the star are not expected to change 
over the pulsation (and the Blazhko) cycle. Hence, if the abun- 
dances are accurately determined at one phase in the pulsation 
cycle, they can be of help to determine (or at least constrain) the 
fundamental parameters at other phases. In this paper we also 
selected the optimal phase for determining the abundances of 
the star. This is the first of a series of planned papers devoted 
to a detailed spectroscopic study of RR Lyr. In forthcoming pa- 
pers we will discuss the spectral variations over the pulsation 
and Blazhko cycle of the star. 

2. Observations and spectral data reduction 

A total of 64 spectra of RR Lyr were obtained between June 
26th and August 27th, 2004 with the Robert G. Tull Coude 
Spectrograph on the 2.7-m telescope of McDonald Observatory. 
This is a cross-dispersed echelle spectrograph yielding a two- 
pixel resolving power R ^ 60000 for the configuration used here. 
Table[T]lists the observing time, exposure time, the S/N ratio per 
resolution element for each acquired spectrum, and the phases 
in the pulsation and Blazhko cycles. For the determination of 
the phases we used the ephemerides derived by Kolenberg et al. 
(2006). To minimize smearing of the spectral features by pul- 
sation, each spectrum was limited to an exposure time of 960 
seconds. Two spectra have shorter exposures as a result of being 
stopped due to cloud. The signal-to-noise ratio (SNR) per res- 
olution element of the obtained spectra varies according to the 
brightness of the star (given the fixed integration time) and the 
weather conditions during the observation. Spectra ID3 19-328 



were inadvertently taken at the wrong blaze angle and thus have 
poor S/N ratios. 

Bias frames and flat-field frames were obtained at the start of 
each night, and Th-Ar spectra were observed frequently during 
each night for calibration purposes. The spectra were reduced 
using the Image Reduction and Analysis Facility (IRAlQ Tody 
1993). Each spectrum, normalised by fitting a low order poly- 
nome to carefully selected continuum points, covers the wave- 
length range 3633-10849 A, with several gaps between the or- 
ders at wavelengths greater than 5880 A. 

The normalisation of the hydrogen lines was a crucial point 
since we used them as source for the derivation of the T e g . Ha 
was not covered by our spectra, and the orders adjacent to H/3 
were affected by a spectrograph defect which hindered a proper 
normalization. We were able to perform a reliable normalisa- 
tion of the Hy line using the artificial flat-fielding technique de- 
scribed by Barklem et al. (2002). This approach assumes that the 
relation between the blaze shapes of the different echelle orders 
is a smoothly changing function of the order number. On this ba- 
sis one can establish the apparent blaze shapes by fitting polyno- 
mials to continuum points in several orders above and below the 
hydrogen line. On the subsequent step a 2-D polynomial surface 
is fitted to these empirical blaze functions, and the continuum in 
the orders containing the Hy line is determined by interpolation. 

This normalisation procedure was performed on two spectra 
taken near the optimal phase. We used a surface fit to 3-4 or- 
ders on both sides of the broad hydrogen line to determine the 
continuum in the Hy spectral orders. The accuracy of this nor- 
malisation technique is attested by a good agreement between 
the normalised overlapping Hy profiles and by the lack of dis- 
crepancy between observations of RR Lyr obtained at similar 
pulsation phases. 

Simultaneously with the spectroscopic campaign, we ob- 
tained photometric data in Johnson V through a multi-site cam- 
paign. The photometric data were published by Kolenberg et al. 
(2006). They were used for accurate determination of the pulsa- 
tion frequencies and phases in the pulsation and Blazhko cycle. 

3. The models 

3.1. The pulsation model 

In order to determine the dynamical properties of the atmosphere 
at the most quiescent phases theoretically, we used the so-called 
Vienna Nonlinear Pulsation Code. For the unperturbed starting 
model we utilized the values of 0.65M Q and 50L Q together with 
r ef f = 6600K and a typical Pop II chemical composition of 
Y - 0.239 and Z = 0.001 which lead to a limit cycle with the 
observed period of pulsation. The kinetic energy of the atmo- 
sphere (defined as the part from the photosphere where r = 2/3 
to the outer boundary of our model, see Fig. [T} shows 2 local 
minima, where the first - roughly at phase 0.35 - corresponds 
to the phase of maximum photospheric radius. The flow ceases 
and the whole envelope starts to contract again to reach its min- 
imal radius shortly after minimum light. Although velocities in 
the atmosphere are low, we are aware that the atmosphere is not 
static at any point during the stars pulsation, which is contrary 
to what is assumed in most model atmosphere codes. But there 
are phases where the dynamical effects are smaller. These are the 



1 IRAF (http://iraf.noao.edu/) is distributed by the National 
Optical Astronomy Observatory, which is operated by the Association 
of Universities for Research in Astronomy (AURA) under cooperative 
agreement with the National Science Foundation. 
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Table 1. Basic data of the observations of RR Lyr. 



pulsation model does not take into account the Blazhko modula- 
tion in the star. 



Spectrum ID 


HJD- 


Pulsational 


Blazhko 


SNR per 


Exposure 


number 


2453000 


phase 


phase 


pixel 


time (s) 


087 


183.7953 


0.173 


0.280 


293 


960.000 


088 


183.8147 


0.207 


0.281 


250 


960.000 


089 


183.8271 


0.229 


0.281 


172 


960.000 


091 


183.8441 


0.260 


0.281 


149 


960.000 


119 


184.7293 


0.820 


0.304 


215 


960.000 


120 


184.7434 


0.846 


0.305 


253 


960.000 


121 


184.7559 


0.868 


0.305 


271 


960.000 


122 


184.7683 


0.890 


0.305 


194 


960.000 


124 


184.7865 


0.922 


0.306 


162 


960.000 


125 


184.7988 


0.943 


0.306 


126 


960.000 


126 


184.8121 


0.967 


0.306 


63 


960.000 


158 


185.7401 


0.604 


0.330 


101 


960.000 


159 


185.7523 


0.626 


0.331 


126 


960.000 


160 


185.7645 


0.647 


0.331 


171 


960.000 


161 


185.7768 


0.669 


0.331 


233 


960.000 


163 


185.7937 


0.698 


0.332 


232 


960.000 


164 


185.8059 


0.720 


0.332 


259 


960.000 


165 


185.8182 


0.741 


0.332 


113 


960.000 


166 


185.8304 


0.763 


0.333 


210 


960.000 


168 


185.8465 


0.792 


0.333 


228 


960.000 


169 


185.8587 


0.814 


0.333 


290 


960.000 


170 


185.8709 


0.834 


0.334 


236 


960.000 


171 


185.8831 


0.856 


0.334 


216 


960.000 


173 


185.8989 


0.885 


0.334 


236 


960.000 


174 


185.9112 


0.905 


0.335 


215 


960.000 


175 


185.9234 


0.928 


0.335 


233 


960.000 


176 


185.9356 


0.948 


0.335 


290 


960.000 


178 


185.9512 


0.977 


0.336 


357 


960.000 


204 


186.7293 


0.349 


0.356 


190 


960.000 


205 


186.7425 


0.372 


0.356 


138 


960.000 


206 


186.7548 


0.394 


0.356 


203 


960.000 


207 


186.7670 


0.416 


0.357 


46 


960.000 


209 


186.7876 


0.452 


0.357 


131 


960.000 


210 


186.8007 


0.475 


0.358 


77 


960.000 


251 


187.7208 


0.098 


0.381 


359 


960.000 


252 


187.7330 


0.120 


0.382 


319 


960.000 


253 


187.7442 


0.141 


0.382 


120 


769.397 


255 


187.7632 


0.173 


0.382 


124 


960.000 


256 


187.7803 


0.203 


0.383 


160 


960.000 


257 


187.7923 


0.226 


0.383 


87 


797.949 


258 


187.8060 


0.249 


0.383 


272 


960.000 


260 


187.8226 


0.278 


0.384 


322 


960.000 


261 


187.8370 


0.303 


0.384 


214 


960.000 


262 


187.8504 


0.327 


0.385 


288 


960.000 


263 


187.8634 


0.349 


0.385 


216 


960.000 


319 


243.6989 


0.854 


1.824 


27 


960.000 


322 


243.7688 


0.976 


1.826 


21 


960.000 


323 


243.7832 


0.002 


1.826 


18 


960.000 


324 


243.7955 


0.024 


1.827 


18 


960.000 


326 


243.8123 


0.053 


1.827 


21 


960.000 


327 


243.8262 


0.078 


1.827 


25 


960.000 


328 


243.8395 


0.101 


1.828 


13 


960.000 


363 


244.6062 


0.454 


1.847 


231 


960.000 


364 


244.6201 


0.478 


1.848 


251 


960.000 


365 


244.6323 


0.500 


1.848 


258 


960.000 


366 


244.6445 


0.521 


1.848 


238 


960.000 


368 


244.6622 


0.553 


1.849 


249 


960.000 


369 


244.6744 


0.574 


1.849 


181 


960.000 


370 


244.6870 


0.596 


1.849 


173 


960.000 


371 


244.6992 


0.618 


1.850 


171 


960.000 


373 


244.7168 


0.648 


1.850 


152 


960.000 


374 


244.7290 


0.670 


1.851 


122 


960.000 


375 


244.7412 


0.693 


1.851 


90 


960.000 


376 


244.7535 


0.713 


1.851 


69 


960.000 



The first column shows the spectrum ID number and the second the 
Heliocentric Julian Date (HJD-2453000) at the beginning of the expo- 
sure. The third and fourth columns show the pulsational and Blazhko 
phase respectively for each observation. The fifth column shows the 
SNR per pixel calculated at ~5000A., while the sixth column shows 
the exposure time, in seconds, for each observation. The spectra num- 
ber 253 and 257 have a lower exposure time, because the observation 
was stopped due to the presence of thick clouds. Spectrum number 260 
is the one analyzed in detail in this paper. 



phases we are interested in for abundance analysis. Note that the 



3.2. The model atmosphere 

To compute model atmospheres of RR Lyr, we employed the 
LLmodels stellar model atmosphere code (Shulyak et al. 2004). 
For all the calculations, Local Thermodynamical Equilibrium 
(LTE) and plane-parallel geometry were assumed. Both these as- 
sumptions need to be discussed in our particular case. 

The LTE assumption may be questionable due to low plasma 
density and shock waves in the RR Lyr atmosphere. To reduce 
the uncertainties of spectroscopic analysis we implement our 
model atmospheres at phases where the dynamical effects in the 
star's atmosphere are expected to be small. We refer the reader to 
the next section for more explanations. Ignoring LTE may also 
lead to systematic errors in the abundance analysis. However, 
a detailed non-LTE analysis is beyond the scope of the present 
paper. 

Atmospheres of giants are extended due to large radii and 
thus sphericity effects may become important for the stellar at- 
mosphere modelling. For instance, based on detailed model at- 
mosphere analysis Heiter & Eriksson (2006) recommended us- 
ing spherically symmetric models for abundance analysis of 
stars having logg < 2 and 4000 K < T eS < 6500 K. Taking 
into account the estimated gravitational acceleration of RR Lyr 
of logg = 2.4 (see next section) one can expect the spheric- 
ity effects to be small enough to not significantly influence 
the abundance analysis. Indeed, the estimated errors of the Fe 
abundance based on theoretical Fei lines presented in Heiter 
& Eriksson (2006) do not exceed 0.1 dex for the model with 
r e fT = 6500, log g = 1 . A log g = 2.4 value for RR Lyr thus justi- 
fies the use of plane-parallel model atmospheres for abundance 
analysis. 

We used the VALD database (Piskunov et al. 1995; Kupka et 
al. 1999, Ryabchikova et al. 1999) as a source of atomic line 
parameters for opacity calculations with the LLmodels code. 
Finally, convection was implemented according to the Canuto 
& Mazzitelli (1991a,b) model of convection (see Heiter et al. 
(2002) for more details). 

4. The optimal phase 

During the pulsation cycle the spectral lines of RR Lyr change 
dramatically. In particular, when a shock wave passes through 
the atmosphere, it is possible to observe line broadening, line 
doubling, line disappearance, and sometimes even line emission 
(Preston et al. 1965, Chadid & Gillet 1996, Chadid et al. 2008, 
Preston 2009). This is obviously the sign of a very non-quiescent 
atmosphere, that in principle cannot be modelled with a static 
model atmosphere, such as ATLAS and LLmodels. In practice a 
model atmosphere code that is able to realistically model the at- 
mosphere of a variable star such as RR Lyr, given both its chem- 
ical and pulsational peculiarity, does not yet exist. Therefore we 
are forced to use a static model atmosphere. In order to anal- 
yse the star in the most consistent way, we decided to study 
RR Lyr when its atmosphere is as close as possible to that of 
a non-variable star. 

In the past, RR Lyrae stars were always analysed through 
spectra obtained close to the phase of minimum light because 
it was believed that this was the phase at which the star's at- 
mosphere is "at its quietest". The adoption of static model at- 
mospheres for the analysis of RR Lyrae stars was then justi- 
fied through showing that the results obtained at different phases 
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were all in agreement with each other. At the "quiescent phase" 
we generally suppose that a) there are no shock-waves or any 
other fast plasma motions in the atmosphere that distort the line 
profiles and b) pressure stratification is as close as possible to its 
hydrostatic analog. 

Picturing a homogeneously pulsating sphere, the most quies- 
cent phase is associated with both phases of minimum and max- 
imum radius. In case of an oversimplified pulsation model, these 
also correspond to the minimum and maximum light in the star's 
luminosity variation. At these extreme positions in the pulsation, 
the atmosphere comes to a halt and, due to zero gas velocity, 
the pressure stratification at these phases is closest to the hydro- 
static case (but not necessarily the same!). Furthermore, since 
the plasma velocities are negligible, the kinetic energy should 
be zero as well. For RR Lyrae stars, and especially those of 
type RRab (fundamental mode pulsators) with strongly nonlin- 
ear light curves (resembling a saw-tooth function), the phase of 
maximum light is short-lived and known to be accompanied by 
shock- waves. Thus, in previous works the quiescent phase was 
associated with the phase of minimum light. However, the as- 
sumption that all atmospheric layers move synchronously as be- 
ing rigidly bound to each other is wrong. Realistic models show 
that at minimum radius, radiation is blocked in deeper layers 
and ready to migrate into the outer layers and accelerate both 
the photosphere and the atmosphere. 

Looking at both pulsational models and a large sample of 
stellar spectra, obtained at different phases, we found that there 
is a quiescent phase very close to the phase of maximum radius, 
where the radial velocity derived from the metallic lines corre- 
sponds to the stellar gamma velocity. On the light curve this cor- 
responds to a phase on the descending branch of the light vari- 
ation. Figure [T] shows the calculated bolometric light curve,the 
photospheric radius variation, and the atmospheric kinetic en- 
ergy as a function of the pulsation phase. The two vertical lines 
correspond to the phases of maximum stellar radius and of min- 
imum light. As mentioned above the expansion of the model en- 
velope is not homologous. We see different waves running in and 
out - some even steepening to shock waves. 

Minimum light is a less fortunate choice to obtain undis- 
torted line profiles, as parts of the atmosphere still move with su- 
personic speed (e.g., Mach 3 in Fig.|2]). Fig.|2]shows radial plots 
covering approximately the outer 270 (of 400 total) radial mesh 
points in the Vienna pulsation model. For the phase of maxi- 
mum radius and the phase of minimum light the gas velocity u 
(in units of the local sound) speed is plotted. At maximum ra- 
dius (left panel in Fig. 2]) the photosphere starts to move towards 
the model's center, while parts of the envelope still move out- 
wards. All velocities are below sound speed. At minimum light 
(right panel in Pig.[2j» - which occurs before minimum radius - 
we see the transition between super- and subsonic inflow, some- 
times called a shock, at the photosphere. This distorts the spectra 
and makes them less suited for our detailed analysis. 

Note that the most quiescent phase, i.e. the phase of maxi- 
mum radius, is actually very short-lived. A spectroscopic obser- 
vation has to be well-timed (within, say, half an hour) in order 
to catch the spectral lines without distortion. Also, the integra- 
tion times cannot exceed a few percent of the pulsation period, in 
order to avoid smearing of the spectra due to pulsation-induced 
Doppler effect. Integration times not much longer than 15 min- 
utes taken within the appropriate (narrow) phase interval are rec- 
ommended. This implies constraints on the obtained signal-to- 
noise ratio, and therefore, this study could only be done with 
+2-m telescopes. At minimum light the star is not at its quietest, 
and the shock wave associated with the bump phase (Gillet & 



Crowe 1988) close to minimum light will also distort the spec- 
tral line profiles in RRab stars. 




0.1 0.4 0.7 

Pulsation Phase [1] 



1.0 



Fig. 1. RR Lyr model bolometric luminosity (upper panel), pho- 
tospheric radius (middle panel), and atmospheric kinetic energy 
(lower panel) as a function of the pulsational phase. As conven- 
tion phase corresponds to the maximum of the luminosity. The 
two full vertical lines correspond to the phases of maximum ra- 
dius (the phase we declare as most quiescent) and of minimum 
light (usually adopted for the spectroscopic analysis). 



Fig.[T] shows clearly that the phase corresponding to mini- 
mum light occurs before a local minimum of the atmospheric 
kinetic energy, while at the phase of maximum radius the stellar 
atmosphere is very close to the other local minimum of kinetic 
energy. 

This picture is also confirmed by the observations. Figure [3] 
shows the comparison of the line profile of RR Lyr in the region 
around 4500 A and the bisector of the Tin line at ~4501 A as 
observed at the phases close to maximum radius and to minimum 
light. The main difference between the two line profiles is given 
by the line broadening, which is an indicator of the atmospheric 
activity: to a quiet phase correspond narrow spectral lines. 

The spectral line broadening as a function of phase is shown 
in Fig. [4] This figure displays the full width at half maximum 
(FWHM) measured for four strong spectral features as a func- 
tion of the pulsational phase. For each of the four lines we ob- 
tained a clear minimum close to the phase of maximum radius 
(between 0.2 and 0.3). This plot also shows the rapid changes 
in the FWHM when the shock waves occur. The observed peaks 
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Fig. 2. Radial plots showing gas velocity u in units of the sound 
speed (the so-called Mach number) shown for the phase of max- 
imum radius (left) and minimum light (right). 



<D 0.7 



o 




Minimum light 
Maximum radius 



wavelength (A) 



velocity (m/s) 



Fig. 3. Left panel: comparison between a part of the RR Lyr 
spectrum as observed close to the phase of maximum radius 
(full line) and close to the phase of minimum light (dashed line). 
Right panel: comparison between the bisectors of the Ti n line at 
~4501 A as observed close to the phase of maximum radius (full 
line) and close to the phase of minimum light (dashed line). The 
dotted line shows the zero velocity, for comparison. The velocity 
scale is in ms~'. 



to Mach 25 in the highest part of the star's atmosphere (see Fig. 3 
in Fokin et al. 1999). 
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Fig. 4. Full width at half maximum (FWHM) in A as a func- 
tion of pulsational phase measured on three Fe n lines and one 
Ti ii line. The minimum is in correspondance with the phase of 
maximum radius. The typical uncertainty on the FWHM is of 
~0.01 A, but it depends also on the SNR of each spectrum. The 
two peaks visible at about phase 0.65 and 0.9 are due to two 
propagating shock waves (Fokin et al. 1999). 



The line bisectors, shown in Figure 3, also display that the 
line asymmetry changes with the pulsational phase. In particular, 
it has a trend similar to the one shown by the FWHM. This point 
will be described in more detailed in the next paper. 

In practice it is not trivial to predict when the star will be 
at maximum radius exactly, due to the presence of the Blazhko 
effect and of other (possibly longer) periodicities involved in the 
pulsation. To find the appropriate spectrum for analysis, we in- 
vestigated all the available spectra obtained close to the phase of 
maximum radius (determined with the simultaneous photome- 
try) and picked the ones showing the minimum FWHM and line 
asymmetry. In the end the two spectra with numbers 258 and 
260 were obtained very close to the phase of maximum radius 
and showed a comparable line broadening and a very small line 
asymmetry. We decided to perform a detailed analysis of spec- 
trum number 260 because of its higher SNR. Our simultaneous 
photometry (Kolenberg et al. 2006) confirmed that this spectrum 
was recorded around pulsation phase <p = 0.28. 



in FWHM near pulsation phase ~0.65 and ~0.9 have been in- 
terpreted as arising from two shocks, a weaker and a stronger 
shock respectively, propagating through the star's atmosphere 
and producing compression of the turbulent gas (Fokin et al. 
1999). Moreover, Fokin & Gillet (1997) and Fokin et al. (1999) 
showed that their RR Lyr models exhibit very strong shocks^jup 



2 Fokin & Gillet (1997) note that the theoretical velocities and the 
shock amplitudes are very sensitive to model parameters. They used pa- 
rameters different from ours: T cff = 7175 K, L = 62L Q , M = 0.578M o , 
X = 0.7 and Y = 0.299. The model generated with the Vienna Nonlinear 
Pulsation Code shows a maximum, during one pulsation cycle, of 4.7 
Mach outward and -3.1 Mach inward, which correspond to gas veloci- 
ties of 35 and -23 km/s. Due to artificial viscosity we can assume that 



5. Fundamental parameters and abundance 
analysis 

In general a fundamental parameter determination begins from 
a derivation of T e g and \ogg from photometric indexes. For 
RR Lyr this operation is not trivial. In Sect. [4] we assumed that 



the star undergoes stronger shocks than our models. The artificial tensor 
viscosity used broadens the shock region and underestimates any possi- 
ble heating phenomena. The position of the outer boundary condition is 
crucial in the maximum velocity, as shock waves steepen up when they 
run outwards towards smaller density. Fokin's model runs till g = 10~ 14 
while the Vienna models we used stop at about 10~ 10 . 
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the atmosphere of RR Lyr can be at best simplified as "static" 
only at the phase of maximum radius. Therefore it is at this phase 
that we made use of the photometric indices and the static model 
atmosphere grids to determine the fundamental parameters. 

As starting point for our analysis we decided to take the 
parameters derived by other authors who analysed spectra of 
RR Lyr obtained at a similar phase. In particular Takeda et 
al. (2006) derived the fundamental parameters spectroscopically 
from high resolution spectra of RR Lyr, one of them obtained not 
far from the phase of maximum radius. These parameters can be 
taken only as starting point because the star was observed at a 
different Blazhko phase^j We will explore the effect of Blazhko 
modulation on the spectra of RR Lyr in our forthcoming work. 

We used the parameters given by Takeda et al. (2006) 
for their spectrum taken at pulsation phase = 0.36 
( r eff =6040+40 K, log g=2.09+0.1 dex) as our starting point. We 
performed an iterative process to improve and test the parame- 
ters as described in the following. In our analysis, every time any 
of the parameters T e ff , log g, i/ m ; c , or abundances changed during 
the iteration process, we calculated a new model by implement- 
ing the most recently determined quantities. We did the same 
with respect to the abundances. While the results of the abun- 
dance analysis depend upon the assumed model atmosphere, the 
atmospheric temperature-pressure structure itself depends upon 
the adopted abundances. We therefore recalculated the model at- 
mosphere every time the abundances were changed, even if the 
other model parameters were kept fixed. This procedure ensured 
that the model structure was consistent with the assumed abun- 
dances. 



5.1. The effective temperature 

We performed the T e g determination by fitting synthetic line pro- 
files, calculated with SYNTH3 (Kochukhov 2007), to the ob- 
served profile of the Hy line, the only hydrogen line for which 
it was possible to make a reliable normalisation. In the tempera- 
ture range expected for RR Lyr, hydrogen lines are very sensitive 
to temperature variations and depend very little on log g varia- 
tions. In particular, this is expected when the stellar T e g is close 
to its minimum. In the case of RR Lyr the use of hydrogen lines 
as T e ff indicators is very important because these lines describe 
the stellar structure more effectively than any other line, being 
formed in a wide region of the stellar atmosphere, and the line 
wings are free from non-LTE effects. The T e g obtained with this 
procedure is T e g — 6125 ± 50 K. Note that these error bars are 
what we obtain from the fitting of the hydrogen line profile. As 
there are model uncertainties that are not taken into account in 
this fitting procedure, these error bars are probably underesti- 
mated. Figure[5]shows the comparison between the observed Hy 
line profile and the synthetic profiles calculated with the adopted 
stellar parameters, as well as the synthetic profiles obtained in- 
creasing and decreasing T e ff by 50 K. 

Figure [5] shows a clear wavelength shift of the observed Hy 
core relative to its wings and to the metal lines formed in deeper 
layers. This is due to a velocity gradient in the outer layers of the 
star (Van Hoof & Struve 1953; Mathias et al. 1995 for RR Lyr). 

Another spectroscopic indicator for T e g is given by the anal- 
ysis of metallic lines. In particular, T e g is determined by elimi- 
nating the correlation between line abundance and line excitation 
potential (y e xcit) for a given ion/element. This procedure can lead 
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Fig. 5. Comparison between the observed Hy line profile (solid 
line) and synthetic profiles calculated with the final adopted 
r e ff= 6125 K (red line) and with r e( f=6175K (green line) and 
T e ff =6075 K (blue line). The red line agrees very well with the 
observed spectrum. 



to erroneous parameters, in particular for stars such as RR Lyr 
where non-LTE effects could be large and the v m [ c determination 
is complicated by the pulsation (see Sect. 



6.1 1. For this reason 



3 Fundamental parameters obtained at different Blazhko phases are 
not necessarily equal to one another, although they are obtained at the 
same pulsation phase. 



we decided not to take this indicator into account in our analy- 
sis, but to use it just as a check of the adopted T e ff value. This 
check pointed towards values comparable to the one we found 
from the Hy line profile fitting. 

5.2. The surface gravity 

The surface gravity can be derived using two independent meth- 
ods: based on line profile fitting of gravity-sensitive metal lines 
with developed wings, and based on the ionisation balance of 
several elements. For RR Lyrae stars the Mg i lines are the most 
suitable lines for the first method. Described in Fuhrmann et al. 
(1997), the first method assumes that the wings of the Mgi lines 
at AA 5167, 5172, and 5183 A are very sensitive to logg varia- 
tions. We decided to use this method only as a check on log g and 
not for parameter determination. The first reason for doing this 
is the large uncertainty in the u sin i and i/ mac ro values that, given 
the available spectral resolution and SNR, could not be precisely 
determined. Another reason is the slight line asymmetry (clearly 
visible for the strong lines, such as Mg i), which makes precise 
line profile fitting impossible. A third reason is the uncertainty 
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in the w m j c value, that will be discussed in detail in Sect. |6.1| To 
make the comparison even more difficult there is the fact that the 
Mg abundance is about 1 dex below solar, making the wings of 
these lines not very pronounced. 

The second method for surface gravity determination uses 
the assumption of ionisation equilibrium, but this method is 
extremely sensitive to the non-LTE effects present for each 
ion/element. Since we are not able to use the line profile fitting 
of the Mg i lines with developed wings, we rely on the ionisa- 
tion equilibrium to determine \ogg, checking the obtained re- 
sult with the Mg i lines. In adopting the ionisation equilibrium, 
for some elements we took into account also the non-LTE cor- 
rections predicted by various authors on some specific elements 
in solar-type metal-poor stars. From the ionisation equilibrium 
we obtained logg=2.4+0.2. This value was derived using only 
the lines with an equivalent width of less than 75 mA, in or- 
der to minimise both the non-LTE effects and the uncertainty 
on the u m \ c (both more pronounced for the strong lines), and to 
assure a large enough number of lines, in particular for iron. We 
checked this value against the observed profiles of the Mg i lines 
with developed wings and obtained a good agreement, in partic- 
ular when we adopt a depth-dependent f m ; c (see Sect. 6.1 1. For 
this comparison we used the line parameters of the Mg i lines 
adopted by Ryabchikova et al. (2009). 

Our value for the surface gravity is supported by the ion- 
isation equilibrium for Fei/Fen and few other elements, such 
as Si i/Si n and Ti i/Ti n. For Ca and V we do not obtain ionisa- 
tion equilibrium, even within the error bars, but we have mea- 
sured just one line for both Ca n and V i. Taking into account the 
non-LTE corrections for Cai (~ +0.1-0.2) and Can (almost in 
LTE) given by Mashonkina et al. (2006), ionisation equilibrium 
is achieved for that element. In the case of chromium, several Cr i 
and Cr n lines have theoretically calculated oscillator strengths, 
which may influence the final abundance results. 

Since RR Lyr's effective temperature is too low to show 
enough He lines (though they are detected - see Preston 2009), 
we are unable to measure the atmospheric He abundance. 
Ryabchikova et al. (2009) tested the effect of a strong He de- 
pletion in the atmosphere of the solar type star HD 49933, con- 
cluding that a depleted He abundance would affect only the log g 
determination and leave T e s unchanged within the error bars of 
0.2 dex. We also tested the effect of a He overabundance for 
RR Lyr. If we assume X = 0.5 and Y = 0.49 we register a gen- 
eral abundance decrease, e.g. Fe decreases 0.3 dex. In addition, 
we observe a variation of the pressure sensitive features such as 
the Mg i lines with extended wings leading to changes in log g 
that do not exceed our error bars. 



5.3. LTE abundance analysis 

Our main source for the atomic parameters of spectral lines is the 
VALD database. LTE abundance analysis was based on equiva- 
lent widths, analysed with a modified version (Tsymbal 1996) 
of the WIDTH9 code (Kurucz 1993). We opted for equivalent 
widths because of the small line asymmetry and of the uncer- 
tainty on the form of the microturbulent velocity, making the 
synthetic line profile fitting more uncertain. We intend to anal- 
yse the other collected spectra of RR Lyr in the same consistent 
way. These spectra show a much more pronounced line asym- 
metry, therefore they will be analysed mostly through equivalent 
widths. 

In total about 700 lines were measured with equivalent 
widths, but after a check against both the solar spectrum and 



the spectrum of HD 49933 (Ryabchikova et al. 2009) we chose 
to keep 617 lines of 26 different elements and 32 different ions. 
We also tried to keep a set of lines uniformly distributed over 
the range of equivalent widths, wavelength, and excitation po- 
tentials, in particular for important ions such as Fe i for which 
we kept 284 lines. We used nearly all unblended spectral lines 
with accurate atomic parameters, except lines in spectral re- 
gions where the continuum normalisation was too uncertain. 
RR Lyr shows a strong underabundance for almost every mea- 
sured ion/element. For this reason, it was possible to give only 
approximate abundance values or upper limits to the abundance 
of some ions for which it was not possible to calculate the equiv- 
alent width because of very weak lines. For these measurements 
we used synthetic line profile fitting, since these lines were too 
shallow to show both any visible line asymmetry and u m ; c de- 
pendence. 

Microturbulence was determined by minimizing the correla- 
tion between equivalent width and abundance for several ions. 
We used mainly Fe i lines since this is the ion that provides the 
largest number of lines within a wide range in equivalent widths, 
but the correlations obtained with Ti i, Ti n, Cr i, Cr n, Fe n, and 
Ni i were also taken into account. Using all the available lines we 
could not find a unique value for the v m j c able to completely re- 
move the correlation between equivalent width and abundance. 
In particular it was possible to remove effectively this correla- 
tion using only the lines with an equivalent width of less than 
75 mA, while the stronger lines exhibited a steep abundance in- 
crease with increasing equivalent width. Using only the lines 
with an equivalent width of less than 75 mA we obtained a v^c 
of 2.4+0.3^8-', while using all available lines we obtained 
a value of 3.1+0.5kms . In this last case the plot shows that 
the correlation is only statistically minimised (see lower plot of 
Fig.|6}. 

In the literature several authors mentioned the possibility of a 
depth-dependent i/ m ; c for RR Lyrae stars (Clementini et al. 1996, 
Takeda et al. 2006). Therefore, given the impossibility of finding 
a clear i> m i c value, we decided to derive the profile of a depth- 
dependent v m \ c using the available measured lines. This part of 
the work will be discussed in Sect. 16.11 

The full set of derived abundances, adopting both a constant 
and a depth dependent u m j c is shown in Table [2] The last column 
of Table [2] gives the solar abundances by Asplund et al. (2005) 
for comparison. 

The stellar metallicity (Z) is defined as follows: 



£ a > 3 m fl 10'°M°-> 



(1) 



where a is the atomic number of an element with atomic mass 
m a . Making use of the abundances obtained from the performed 
analysis and assuming the depth-dependent v m [ c , we derived a 
metallicity of Z = 0.003 ± 0.002 dex. For elements that were 
not analysed we adopted solar abundances from Asplund et al. 
(2005). However, if we assume an underabundance of -1.0 dex 
for all elements that were not analysed, excluding H and He, the 
resulting Z value remains practically unchanged. Note that in the 
pulsation model we used a metallicity value of Z = 0.001. This 
is on the lower border of the value we obtain from our spec- 
trum. However, we have to keep in mind that we do not take 
into account NLTE effects in our analysis, and that these ef- 
fects may influence our result. In general, for metal-poor stars 
such as RR Lyr, the NLTE correction is negative, meaning that 
the abundance in NLTE is lower, so the real Z would be lower 
too (from 0.003 more towards 0.001). The effect of NLTE on 



8 



K. Kolenberg et al.: An In-Depth Spectroscopic Analysis of the Blazhko Star RR Lyr 



Table 2. LTE atmospheric abundances for RR Lyr. 



Ion 




RR Lyr 






Sun 




logWJVto,) 


log(JV/JV M ) 


n 


Remarks 


log(A7iV tol ) 


Ci 


-4.84±0.()4 


-4.78±0.07 


4 




-3.65 


Ni 


-4.77±0. 1 1 


-4.77±0.11 


2 


S 


-4.26 


Oi 


-3.98±0.04 


-3.97±0.04 


6 


S 


-3.38 


Nai 


-7.38±0.05 


-7.38±0.05 


2 




-5.87 


Mgi 


-5.54±0.10 


-5.48±0.14 


8 




-4.51 


Ah 


-6.59±0.26 


-6.57±0.26 


9 


S 


-5.67 


Si i 


-5.56±0.11 


-5.55±0.11 


10 




-4.53 


Sin 


-5.61±0.04 


-5.54±0.08 


3 




-4.53 


Si 


-5.94±0.08 


-5.93±0.08 


7 


S 


-4.90 


Ki 


-7.81 


-7.67 


1 




-6.96 


Cai 


-6.98±0.05 


-6.90±0.09 


24 




-5.73 


Ca ii 


-6.71 


-6.69 


1 




-5.73 


Sc ii 


-10.17±0.13 


-10.18±0.09 


17 




-8.99 


Til 


-8.33±0.08 


-8.32±0.07 


22 




-7.14 


Tin 


-8.21±0.19 


-8.22±0.13 


56 




-7.14 


Vi 


-9.64 


-9.62 


1 




-8.04 


Vn 


-9.22±0.09 


-9.19±0.09 


8 




-8.04 


Cri 


-7.98±0.09 


-7.94±0.09 


18 




-6.40 


Cm 


-7.70±0.11 


-7.66±0.10 


27 




-6.40 


Mm 


-8.47±0.15 


-8.42±0.12 


7 




-6.65 


Fei 


-6.07±0.12 


-6.03±0.11 


284 




-4.59 


Fe n 


-5.93±0.13 


-5.89±0.10 


47 




-4.59 


Co i 


-8.22 


-8.21 


1 




-7.12 


Nil 


-7.35±0.08 


-7.33±0.08 


38 




-5.81 


Cui 


-9.74 


-9.73 


1 




-7.83 


Zn i 


-9.01±0.01 


-8.99±0.01 


2 




-7.44 


Gai 


-10.50 


-10.50 


2 


UL/S 


-9.16 


Rbi 


-9.75 


-9.75 


1 


UL/S 


-9.44 


Sri 


-10.47 


-10.46 


1 




-9.12 


Srn 


-10.50 


-10.48 


1 




-9.12 


Yn 


-11.31±0.09 


-11.30±0.09 


10 




-9.83 


Zrn 


-10.56±0.07 


-10.53±0.06 


4 




-9.45 


Nb n 


-11.50 


-11.50 


1 


UL/S 


-10.62 


Mo i 


-10.50 


-10.50 


1 


UL/S 


-10.12 


Pdi 


-11.00 


-11.00 


1 


UL/S 


-10.35 


Ban 


-11.23 ±0.22 


-11.35±0.08 


5 




-9.87 


La ii 


-12.05±0.02 


-12.05±0.02 


3 




-10.91 


Ce n 


-11.74±0.08 


-11.73±0.08 


8 




-10.46 


Pru 


-12.36±0.10 


-12.35±0.10 


16 


S 


-11.33 


Nd n 


-11.92 


-11.91 


1 




-10.59 


Smii 


-12.00±0.07 


-12.00±0.07 


10 


S 


-11.03 


Eun 


-12.57 


-12.55 


1 




-11.52 


Gdn 


-12.00±0.07 


-12.00±0.07 


5 


S 


-10.92 


Tbn 


-12.25±0.35 


-12.25±0.35 


2 


S 


-11.76 


Dy n 


-11.82±0.19 


-11.80±0.19 


6 


S 


-10.90 


Hon 


-12.75±0.35 


-12.75±0.35 


2 


s 


-11.53 


Ern 


-11.61 


-11.60 


1 




-11.11 


Tmii 


-12.50 


-12.50 


1 


UL/S 


-12.04 


Lun 


-12.83 ±0.15 


-12.83±0.15 


3 


S 


-11.98 


Hfn 


-11.50 


-11.50 


1 


UL/S 


-11.16 


Pbi 


-11.00 


-11.00 


1 


UL/S 


-10.04 


Thu 


-12.90 


-12.90 


1 


UL/S 


-11.95 


Trf 




6125 K 






5777 K 


logs 




2.40 






4.44 




3.1 km s" 1 


Depth dep. 






0.875 



Error bar are estimates based on the internal scatter from the number 
of analysed lines, n. The third column gives the atmospheric abun- 
dances in case of a polynomial depth dependent microturbulence ve- 
locity. The last column gives the solar abundance values from |Asplund| 
et al. (2005 ). The column indicated as "Remarks" shows whether the 
given abundance value is an upper limit (UL) and/or was obtained with 
synthetic spectra (S). 



the abundances was illustrated by, e.g., Takeda et al. (2006) for 
oxygen. Using the solar Fe abundance value of Asplund et al. 
(2005), we obtain [Fe/H]=-1.41 + 0.11 for RR Lyr. This is in 
good agreement with the value [Fe/H]=-1.39 obtained by Beers 
et al. (2000), who list typical errors of 0.1-0.2 dex. 

We used spectral synthesis with SYNTH3 to check the hy- 
perfine structure (hfs) effects on the abundance determination of 
Mn, Cu, Zn, Ba, and Pr. For each measured line of each of these 
elements hfs effects are less than 0.01 dex, except for the Ban 
line at A 6141 A for which the hfs correction is -0.1 dex, bring- 



ing the line abundance closer to the mean Ba abundance. The 
hfs calculations for barium were taken from Mc William (1998), 
who does not list the parameters for the Ba line at A 6496 A, but 
for this specific line we do not expect any significant hfs effect 
(Mashonkina & Zhai 2006). 

The abundance uncertainties given in Table [2] are the stan- 
dard deviation from the mean abundance (hence no uncertain- 
ties are given if the abundances were derived from a single line). 
More realistic error bars for each element/ion can be found in 
Ryabchikova et al. (2009) where a rigorous derivation of the 
abundance uncertainties is given on the basis of the adopted un- 
certainties on the stellar parameters. This direct comparison is 
possible because RR Lyr and HD 49933 have a similar T e ff and 
in particular similar values of the uncertainties on both T e s and 

Given the quality of the data and the slight line asymmetry it 
was not possible to give definite values for both u sin i and f ma cro, 
but just to constrain their values. We obtained that v sin i lies be- 
tween and 9kms~ I , and v m . dcm between 6 and llkms -1 . In 
particular with the minimum given v sin i we obtained the max- 
imum f maclo and vice versa, as both effects contribute to line 
broadening. Note that the spectral lines cannot be fit with only 
the effects of rotational broadening (v sin z). The constraints on 
v sin i that we obtained are in good agreement with the one ob- 
tained by Kolenberg (2002) through analysis of the line profile 
variations of the star. They also are in accord with the findings by 
Peterson et al. (1996), who measured the line widths for 27 RR 
Lyrae (of which 8 RRab) variables via cross-correlation analy- 
sis. They estimated an upper limit of 10 km/s for vsini in all 
cases. 



6. Discussion 

6.1. A depth-dependent microturbulent velocity 

As previously mentioned, we calculated the profile of a depth- 
dependent v m i c on the basis of the available equivalent widths 
obtained for several Fe i lines spanning a large range of values. 
A depth-dependent v m i c was previously suggested by several au- 
thors, such as Takeda et al. (2006), who found evidence that 
"strongly suggest that the microturbulence increases with height 
in the atmosphere of RR Lyrae stars, and that a simple applica- 
tion of the v m i c value derived from weak/medium-strength lines 
to stronger lines may result in an overestimation of the abun- 
dances". 

Figure [6] displays the line abundance as a function of the 
measured equivalent width for all the measured Fe i lines in the 
spectrum of RR Lyr, calculated by assuming constant (bottom) 
and depth-dependent (middle) v m \ c . It is clear that the use of a 
constant v mlc leads to an underestimation of the abundance of 
the medium-strength lines and an overestimation for the strong 
lines. In Fig. [6] we included as comparison the line abundance 
as a function of the measured equivalent width in HD 49933 for 
the set of common Fe i lines. This demonstrates that the observed 
behavior does not depend on the set of adopted lines. 

We believe that this phenomenon could be explained both 
by strong non-LTE effects and by a depth-dependent v m i c . 
Generally speaking, non-LTE effects are stronger for deep lines 
compared to shallow lines and adopting line formation in LTE 
would lead to a higher abundance, in agreement with what 
we observe here. Gehren et al. (2001) showed that for solar- 
type stars there is a substantial Fe i underpopulation leading to 
stronger Fe line wings when assumed in LTE, while Fe n is in 
LTE, but Gehren et al. (2001) adopted a model atom for Fe that 
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o constant Vraic = 3.1 km/s 

□ polynomial Vmic 

* HD 49933: Vmic = 1.6 km/s 




equivalent width (mA) 

Fig. 6. Plots of individual abundances for 284 Fe i lines versus 
the measured equivalent width for RR Lyr, adopting a constant 
(open circles) and depth-dependent (open squares) u m ; c . The 
same dependence for 108 common Fei lines in the spectrum of 
HD 49933 is shown by crosses. We applied an arbitrary vertical 
shift for visualisation purposes. 



did not include high-excitation levels. Mashonkina et al. (2009) 
analysed the Fe i/Fe n ionisation equilibrium in four solar-type 
stars and in the Sun concluding that the inclusion of the high- 
excitation levels in the Fe i model atom substantially reduced the 
non-LTE effects. RR Lyr is a metal poor giant for which non- 
LTE effects are expected to be stronger than in solar-type stars. 
If the deviation we register is due to non-LTE effects, it should 
only produce a deviation of about 0.5 dex for the stronger lines. 
We believe that the non-LTE effects are only partially responsi- 
ble for the obtained deviation, because we observe it in all mea- 
sured ions with a similar magnitude and always in the same di- 
rection. These ions include those for which non-LTE effects are 
supposed to be weak, such as Fe n, that shows deviations even 
stronger than those registered for Fe r. Note that non-LTE effects 
work differently for different ions, leading to deviations in both 
directions and with a wide range of magnitudes. For this rea- 
son we believe that a depth-dependent y m j c is mostly responsible 
for the observed deviations. A depth-dependent i/ m ; c is also sup- 
ported by modelling of RR Lyr stars as shown by Fokin et al. 
(1999) who observed pulsation-dependent variations in the mi- 
croturbulent velocity (see their Fig. 3). 

The depth-dependent v m i c profile was obtained fitting the 
line abundance in the plane equivalent width versus line abun- 
dance. We performed this procedure for Fei, given the large 
number of measured lines and then tested the solution with the 
other ions. For each line the code searches the best individual 
line abundance, assuming a certain dependence of the micro- 
turbulent velocity on the atmospheric depth. This dependence 
is varied to minimise the dispersion between the observed and 
the theoretical line widths over the whole set of measured spec- 
tral lines. As the equivalent widths come from lines which are 
formed in a small fraction of the stellar atmosphere, it is impos- 



4 The code adopts an LTE plane-parallel model atmosphere cal- 
culated using subroutines of the SynthV code (Tsymbal 1996) and 
the DUNLSF minimisation procedure of the IMSL numerical libraries 
package. 



sible to obtain a u m i c value at each atmospheric depth. For this 
reason it was necessary to speculate about the analytic form of 
the i/ m i c dependence on the atmospheric depth. We tested both 
a step-like function and a low-degree polynomial function. We 
chose the latter due to the unrealistically steep and large step 
needed for the step-like function. 

Figure [7] shows the profile obtained for the depth-dependent 
v m \ c in comparison with the sound speed calculated by 
LLmodels. According to the results of Fokin et al. (1999) the 
v mic should always be subsonic, due to the presence of strong 
dissipation effects. In Fig. [7] this is the case, except for the re- 
gion between logrsooo equal to -3.5 and -5.5 v m { c where i/ m i c 
becomes supersonic. This, however, should be interpreted with 
caution since the numeric calculation of the sound speed with the 
LLmodels code may suffer from accuracy loss in uppermost lay- 
ers where the thermodynamic variables (like pressure and den- 
sity) are slowly changing functions with atmospheric depth. Also 
notice that our empirical estimation of the u m ; c may contain sys- 
tematic uncertainties, and thus the supersonic regime shown in 
Fig.|7]may have little to do with reality. Nevertheless, the general 
behavior of i/ m ; c with depth (i.e., strong increase in superficial 
layers) plauseably reflects a real physical phenomenon. 



microturbulence velocity 

sound speed 



25 - 

t/3 

^20 - 




L °g CW 

Fig. 7. Comparison between the profile of the microturbulent ve- 
locity (full line) and of the sound speed (dashed line) as a func- 
tion of atmospheric depth. 



Assuming a depth-dependent i/ m ; c the abundances of most 
ions change, in particular the ones with a large number of strong 
lines. As expected, for most of the ions the abundance increases 
instead of decreasing. This phenomenon is explained by the fact 
that our adopted value of constant i/ m ; c is overestimated because 
it was obtained by minimising statistically the correlation abun- 
dance versus equivalent width. Because we have a large num- 
ber of middle-strength lines, a small number of strong lines 
and an overestimated y m ; c we obtained an underestimated abun- 
dance (see the bottom plot of Fig. [6]). When applying the depth- 
dependent v m i c , the middle-strength lines (responsible for the 
underestimated abundance) yield an abundance similar to the 
shallow and strong lines, thus leading to a higher final abun- 
dance. This effect is visible for ions with a large number of lines, 
such as Fe i. Instead, for elements with few lines and some of 
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them very strong, such as Ban, the final abundance decreases, 
with a considerable decrease of the standard deviation as well. 

6.2. Comparison with previous determinations 

Since we propose a different approach to the analysis of RR 
Lyrae stars with the aim of determining their physical param- 
eters, it is important to compare our results with what was pre- 
viously obtained by other authors. The main publications on this 
topic containing an abundance analysis and parameter determi- 
nation of RR Lyr itself are Clementini et al. (1995), Lambert et 
al. (1996), and Takeda et al. (2006). 

Clementini et al. (1995) analysed several field RR Lyrae 
stars to derive broad band photometric calibrations for both fun- 
damental parameters and metallicity. Their spectra of RR Lyr 
covered a large wavelength range with a moderate resolution 
(R-38 000) and a high SNR (probably assumed per pixel) of 
about 460. They observed RR Lyr at pulsational phases of 0.70 
and 0.73, around minimum light. Since they did not find any sys- 
tematic difference between the two spectra they summed them 
to get one high-SNR spectrum. They compare the r er y value ob- 
tained from broad band photometry by several other authors be- 
tween 1975 and 1994, and, taking into account these compar- 
isons, they build their own calibration leading to a T e ff value 
of 6222+115 K. They derive the surface gravity from the stel- 
lar mass and radius, obtaining a logg value of 2.8±0.2dex. The 
microturbulent velocity was derived in the usual way (minimi- 
sation of the correlation between line abundance and equivalent 
widths) obtaining a value of 4.2+0.2 km s . Clementini et al. 
(1995) also mention the possiblity of a depth-dependent u m i c , but 
they concluded that, if present, a u m ; c depth-dependency is only 
very small and affects the abundances by just 0.1 dex (0.2 dex 
for the elements with a large number of strong lines). We con- 
firm this estimate with our present work (see Table|2]i. Most of 
their abundances are in LTE, exept for oxygen and sodium that 
were analysed in non-LTE. Table[3]shows a comparison between 
the abundances obtained by Clementini et al. (1995) and the ones 
derived in this work. The comparison shows a rather good gen- 
eral agreement between the two sets of abundances, where for 
only Al i, Sc a, Mn i and Zn i we register a small disagreement. 

Lambert et al. (1996) observed a set of RR Lyrae stars with 
a spectral resolution of 23 000 to obtain narrow and broad band 
photometric calibrations for fundamental parameters and metal- 
licity. RR Lyr was observed at eight different phases, one of 
them close to the phase of minimum light and one close to 
the phase of maximum radius. The spectrum obtained close to 
the phase of maximum radius has a SNR (assumed per pixel) 
of ~100. For each phase they derived the fundamental parame- 
ters both from photometry (adopting previously existing calibra- 
tions) and from spectroscopy. From the photometric calibration 
they obtained, for the spectrum at phase close to maximum ra- 
dius, T e ff =6350+200 K and log g=2.6±0.2 dex, while from spec- 
troscopy they derived T e g =6200+200 K and logg=2.3±0.2dex. 
Lambert et al. (1996) found a constant u m - ic along the stellar at- 
mosphere, but the most surprising result is that they obtained 
also a rather constant v^c along the pulsational phase (between 
3.6+0.3 km s and 4.4+0.5 km s -1 ). This result is surprising 
given the fact that several pulsation models of RR Lyrae stars 
show large variations of the u m ; c along the pulsational cycle. We 
will study this issue in detail in the forthcoming work. In gen- 
eral the high value obtained by Lambert et al. (1996) shows the 
turbulent motions present in the atmosphere and it fits with the 
large value we also obtained when a constant Vmc is assumed. 
Lambert et al. (1996) derived the Ca and Fe abundance both in 



Table 3. Comparison between the atmospheric ion abundances 
relative to the Sun obtained by Clementini et al. (1995) and in 
this work. 



Ion 




[AWiVH 


s 






This work 




Clementini et al.||1995| 


Oi 


-0.60±0.04 


-0.59±0.04 


-0.69 


Nai 


-1.51±0.05 


-1.51±0.05 


-1.39 


Mgi 


-1.03±0.10 


-0.97±0.14 


-1.08±0.07 


Ah 


-0.92±0.26 


-0.90±0.26 


-1.89 


Si i 


-1.03±0.11 


-1.02*0.11 


-0.92±0.02 


Si ii 


-1.08±0.04 


-1.01 ±0.08 


-1.14 


Cai 


-1.25±0.05 


-1.17±0.09 


-1.07±0.05 


Sen 


-1.18±0.13 


-1.19±0.09 


-1.35±0.12 


Til 


-1.19±0.08 


-1.18±0.07 


-1.20±0.06 


Tin 


-1.07±0.19 


-1.08±0.13 


-1.05±0.18 


Cn 


-1.58±0.09 


-1.54±0.09 


-1.41±0.18 


Crn 


-1.30±0.11 


-1.26±0.10 


-1.33±0.14 


Mm 


-1.82±0.15 


-1.77±0.12 


-1.99±0.12 


Fei 


-1.48±0.12 


-1.44±0.11 


-1.39±0.13 


Fen 


-1.34±0.13 


-1.30±0.10 


-1.39±0.13 


Nil 


-1.54±0.08 


-1.52±0.08 


-1.47 


Zn i 


-1.57±0.01 


-1.55±0.01 


-1.35 


Yn 


-1.48±0.09 


-1.47±0.09 


-1.50 


Ban 


-1.36±0.22 


-1.48±0.08 


-1.39±0.12 


Willie 


3.1kms-' 


Depth dep. 


4.2 km s" 1 



The second and third column show the abundances relative to the Sun 
obtained in this work, while the fourth column lists the abundances rel- 
ative to the Sun published by Clementini et al. (1995). 



LTE and non-LTE, adopting only the lines that are not deeper 
than 100 mA. In LTE they obtained log(/V Fe //V to t)=-6.06 dex and 
log(A r c a /Mot)=-6.89 dex. In non-LTE they obtained a correction 
of about 0.2 dex for Fe i and no correction for Fe n, but, as men- 
tioned in Sect. |6.1| we believe that the non-LTE correction for 
Fe i is much smaller. The non-LTE correction they obtained for 
Ca i was ~0.05 dex. Both the Fe and Ca abundance presented 
by Lambert et al. (1996) are in good agreement with what we 
obtained. There is the possibility that the higher T e g and v m \ c 
relative to what we adopted compensate each other leading to 
values very close to the ones obtained with a lower and u m ; c . 

Takeda et al. (2006) analysed five spectra of RR Lyr with the 
aim to derive spectroscopically the fundamental parameters and 
the abundances of O, Si, and Fe. The spectra were obtained in 
high resolution (R -60 000) with a rather high SNR (-350-400), 
but in a very limited wavelength range intended to cover mainly 
the oxygen triplet at AA -7770 A. One of these five spectra was 
obtained at a pulsational phase close to the phase of maximum 
radius. The fundamental parameters were derived minimising 
the correlation between line abundance and excitation potential 
(r er r) and with the ionisation equilibrium (logg). The hydrogen 
line wings were used as check of the parameters determined, the 
inverse of the strategy we applied. For this phase they obtained 
r eff =6040+40 K, logg=2.09±0.1dex, and y mic =3.0+0.1kms" 1 . 
The v m \ c value obtained by Takeda et al. (2006) is in very good 
agreement with what we have found using a constant y m ; c . They 
also strongly suggest the presence of a depth-dependent v m i c for 
RR Lyrae stars. Takeda et al. (2006) analysed two sets of oxygen 
lines in non-LTE obtaining two values for the O abundance of - 
4.11 dex and -4.00 dex. We believe that the difference between 
their two obtained values is due to the fact that they adopted a 
constant w m j c . The strong infrared oxygen triplet is very sensitive 
to the adopted v m \ c . The O abundance we obtained is very close 
to -4.00 dex assuming both a constant and a depth-dependent 
v m { c . Takeda et al. (2006) also derived the Si and Fe abundance in 
LTE, obtaining respectively -5.93 dex and -5.82 dex. These two 
values do not match very well our results (we obtain a higher Si 



K. Kolenberg et al.: An In-Depth Spectroscopic Analysis of the Blazhko Star RR Lyr 



11 



abundance and a lower Fe abundance). We are not able to explain 
this difference since it cannot be due to the small differences in 
the adopted stellar parameters. 

The effective temperatures determined in these studies may 
not be directly comparable because the analysed spectra were 
obtained at different Blazhko phases A more reasonable compar- 
ison can be done for the surface gravity because it is supposed 
to change less during the pulsation cycle. 

What clearly emerges is a systematic difference between the 
log g obtained from stellar mass and radius and from the ioni- 
sation equilibrium. The ionisation equilibrium leads to a lower 
log g compared to what was derived by the assumed stellar mass 
and radius. Our \ogg value lies in between and it is well vali- 
dated from the fit of the magnesium lines with developed wings. 
This method is independent of both the ionisation equilibrium 
and the assumed stellar mass and radius. 

The other parameters that can be directly compared are the 
elemental abundances. We do not expect those to vary within 
the pulsation cycle, unless the pulsation is able to bring material 
from the inner core up to the stellar atmosphere, a possibility 
that we find unlikely. Comparing the abundances obtained in the 
three previous studies and in our work, we note good agreement 
except for a few exceptions. This means that there is the actual 
possibility that if static models can be applied to RR Lyrae stars, 
they can be applied to many phases (not all!), but we will tackle 
this issue in our next work on the analysis of the other acquired 
spectra of RR Lyr. 

From the point of view of the atmospheric modelling, the 
main difference between earlier work and ours is the adoption 
of an abundance-dependent model atmosphere for an RR Lyrae 
type star. The impact of depth-dependent microturbulent veloc- 
ity and individual chemistry is shown in Fig. [8] As expected, 
strong underabundance makes the temperature in the surface lay- 
ers higher compared to solar or scaled solar models. Introduction 
of the depth-dependent v m \ c affects the strongest spectral fea- 
tures, leading to weaker absorption at their line center. This, in 
turn, decreases the line absorption coefficient and thus leads to a 
further increase of the temperature. For instance, the temperature 
difference between models computed with constant and depth- 
dependent microturbulent velocities (see Fig. 7) can amount to 
as much as 300 K. On the other hand, only the strongest lines 
feel this modification in model temperature structure and thus 
the statistical results of abundance analysis are not affected that 
much, as shown in Table 2. 

In contrast, line profile analysis of strong lines must be 
performed with appropriate models that account for the ef- 
fect of depth-dependent microturbulence. Finally, we do not 
find any noticeable difference in temperature structure between 
the scaled-solar model ([M/H]=-1.5), which corresponds to the 
mean underabundance we found for Fe, and individual abun- 
dance models. This allows us to use scaled abundance models 
to mimic the temperature structure of the star as a first guess. 
Note however that metallicity is the result of abundance analy- 
sis, which is not known in the beginning. This methodological 
difficulty is automatically removed using an iterative procedure 
of abundance analysis as applied in this study. 

7. Conclusions 

Fundamental-mode RR Lyrae (RRab) stars pulsate with high ve- 
locities which, in certain phases, leads to a distortion of their 
spectral features. The aim of the study presented in this paper 
was to determine a set of self-consistent and accurate param- 
eters and abundances for RR Lyr, the prototype and brightest 
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Fig. 8. Comparison of temperature structures of models com- 
puted under different assumptions about microturbulent veloc- 
ity fmic and abundances (see the legend on the plot). Depth- 
dependent v m [ c was taken from Fig. [7] For all the models we 
took r eff = 6125, log g= 2.4. 



member of its class. For this, we had a set of high-resolution 
spectroscopic data at our disposal, obtained with the Robert G. 
Tull Coude Spectrograph on the 2.7-m telescope of McDonald 
Observatory. In order to derive the most reliable abundances for 
the star, we determined the phase in RR Lyr's pulsation cycle 
at which the atmosphere is "at its most quiescent", the phase 
of maximum radius. Our assumptions were strengthened by the 
Vienna Nonlinear Pulsation Code. With observations taken at 
the phase of maximum radius, the fundamental parameters and 
the element abundances of RR Lyr were determined through an 
iterative process. For the determination of the effective temper- 
ature, we used synthetic line profile fitting to the Hy line and 
obtained T e ff = 6125 ± 50 K (error bars from the fitting proce- 
dure only). For the \ogg determination we used the condition of 
ionisation equilibrium and obtained a logg value of 2.4 + 0.2 
at the phase of maximum radius. The LTE abundance analysis 
based on element abundances versus equivalent widths could not 
fit the element abundances with a single microturbulent veloc- 
ity (v m ; c ) value. A depth-dependent v m i c is physically plausible 
for RR Lyrae atmospheres and was previously suggested (but 
never derived or quantified) by several authors. In this work, 
we derived, for the first time, the depth-dependent v m [ c profile 
and quantified the expected abundance variation. With the depth- 
dependent v m \ c we obtained a better agreement between the ele- 
ment abundances and equivalent widths. In general, the adoption 
of a fixed i/ m j c value (which is too high) leads to an underestima- 
tion of the element abundances by 0.06 dex (maximum value). 
Nevertheless, RR Lyr is shown to be underabundant in all heavy 
elements, in agreement with previous studies. 

RR Lyr, the eponym of its class, is one of the best studied 
RR Lyrae stars. However, many intricacies of its pulsation re- 
main poorly understood, and in order to accurately model the 
star, we have to take into account complex physics that we are 
only beginning to uncover. The star is an asteroseismic target 
of the Kepler Mission (Borucki et al. 2009) through the Kepler 
Asteroseismic Science Consortium (KASC, see also Gilliland et 
al. 2010). Besides the fact that the star displays amplitude and 
phase modulation (the so-called Blazhko effect), Kolenberg et 
al. (2010) recently detected, for the first time, the occurrence of 
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half-integer frequencies in the star, i.e., peaks at j/o, |/o, |/o 
etc., with /o the main pulsation frequency . This phenomenon 
("period doubling") may be caused by instabilities in the star 
and connected to the mysterious Blazhko effect (see Szabo et 
al. 2010, in preparation). In order to further explore theoretical 
models of RR Lyr, it is necessary to know the physical parame- 
ters of the star with the highest accuracy possible. 

For this reason, an in-depth analysis of the star's atmospheric 
motions is very timely. We intend to expand this analysis in our 
forthcoming publications of the star. Our study clearly illustrates 
that it is crucial to use the appropriate models to correctly inter- 
pret the spectral data of RR Lyrae stars, and that high-quality 
observations can contribute to improving those models. 
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